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Abstract. The low temperature acoustic and thermal properties of amorphous, 
glassy materials are remarkably similar. All these properties are described 
theoretically with reasonable quantitative accuracy by assuming that the 
amorphous solid contains dynamical defects that can be described at low 
temperatures as an ensemble of two- level systems (TLS), but the deep nature 
of these TLSs is not clarified yet. Moreover, glassy properties were found also in 
disordered crystals, quasicrystals, and even perfect crystals with a large number 
of atoms in the unit cell. In crystals, the glassy properties are not universal, like 
in amorphous materials, and also exhibit anisotropy. Recently it was proposed 
a model for the interaction of two-level systems with arbitrary strain fields 
(Phys. Rev. B 75, 64202, 2007), which was used to calculate the thermal 
properties of nanoscopic membranes at low temperatures. The model is also 
suitable for the description of anisotropic crystals. We describe here the results 
of the calculation of anisotropic glass-like properties in crystals of various lattice 
symmetries, emphasizing the tetragonal symmetry. 



1. Introduction 

The low temperature acoustic and thermal properties of amorphous, glassy materials 
are remarkably similar and they can be explained-to a large extent-by assuming that 
the material contains a large number of dynamic defects. These dynamic defects are 
tunneling systems (TS) which are modeled by an ensemble of two-level systems (TLS) 
0; 0]. Crystals with defects-with a large enough amount of disorder-exhibit also 
glass-like properties, but these properties are not so universal and, even more, they 
are not isotropic-for example the absorption and the propagation velocity change of 
elastic waves depend on the propagation and polarization directions of the wave [3(] . 

A microscopic model of the tunneling systems is still missing. For this reason, 
the study of disordered crystals is especially interesting because it offers additional 
information: in some materials we know quite well which are the entities that tunnel 
between different equilibrium positions. Beside this, the anisotropy of the TLS-sound 
wave interaction in crystals represents another challenge to the interaction models of 
TLSs which requires clarification. 

A free TLS is described in a two-dimensional Hilbert space by the Hamiltonian 




(1) 
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where o~ z and o~ x are Pauli matrices; the parameters A and A are called the asymmetry 
of the potential and tunnel splitting, respectively. The excitation rate of the TLS, e, 
is obtained by diagonalizing ifjLS and this gives e = \/A 2 + A 2 . 

The perturbation of the free TLS, induced by a strain field is H\ = (S/2)a z , where 
in general S = 2jij Sij , with Sij being the strain tensor and 7^ a symmetric tensor that 
characterizes the TLS and its interaction with the strain field-here, as everywhere in 
this paper, we assume summation over repeated indices. In bulk isotropic materials, 
like in amorphous solids, the elastic waves may be decomposed into longitudinally and 
transversely plane waves and the interaction term 5 may be written as 

S = 2 la S a , (2) 

where a is the polarization, I (longitudinal) or t (transversal), of the perturbing wave. 
In this case the absorption rate of a phonon of polarization a and wavevector k, by a 
TLS is 

2tt 2 hk A 2 

rk, CT = -rloTT} jn k ^d(hc a k- e), (3) 

n 2V pc a e z 

where by n^ a , p, c CT , and V we denote the population of the phonon mode, the 
density of the solid, the sound velocity and the volume of the solid, respectively. If 
the number of TLSs of asymmetry and tunnel splitting in the intervals dA and dA 
is P(A,A)dAdA = VP /AdAdA, then integrating T k(T over the TLS ensemble and 
using 7ik,o- = [exp(f3hca-k) — 1]~ , we obtain the standard total phonon absorption 
rate, 

(«) _1 =^-tanh(^), (4) 

where we introduced the short-hand notation K = Tm^ a k / '(V pc a ) . 

Although widely used, the simplified version of the tunneling model exposed above 
cannot account for the anisotropy of the glass-like properties observed in disordered 
crystals. To explain this anisotropy, in Refs. 0; @| we applied the model introduced 
in 0. In this model, a direction t is associated to each TLS [t = (t x ,t y ,t z Y , where 
the superscript ' denotes the transpose of a matrix or vector] and S was written in the 
form [1 H 1 

S = 2T* ■ [R] ■ S, (5) 

with T = [t^.,ty,t^.,2tyt z ,2t z t x ,2t x ty) t and S is the strain tensor in abbreviated 
subscript notations. The 6x6 symmetric matrix [R] contains the TLS-strain field 
coupling constants ry and its structure is determined by the symmetries of the lattice 
in which the TLS is embedded. In Refs. 0; @| we applied this model to calculate the 
glassy properties of disordered cubic (Ca stabilized zirconium 0), trigonal (neutron 
irradiated quartz 10]) and hexagonal crystals. Here we shall report results on 
tetragonal lattices. 



2. Phonon scattering rates in tetragonal lattices 

In the new model, the phonon absorption rate by a TLS is 

r kCT (t) = y ^^|T 4 • [R] ■ S ka \ 2 S(e - M. (6) 

We see that the main characteristic of the TLS-phonon interaction is contained in the 
quantity M^ a (t) = T* • [R] ■ Sko-, which bears an intrinsic anisotropy through the 
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matrix [R] , on which the symmetries of the lattice are imposed [1; 0; @; 0| To calculate 
the average scattering rate of a phonon by the ensemble of TLSs, now we have to 
average not only over the distribution of A and A, but also over the distribution of t. 
In this way we get the total phonon absorption rate, 

| 2irP tanh ( 



2fc B T 



-nk, 



(7) 



To reduce the number of degrees of freedom of the problem, in what follows we shall 
assume that t is isotropically oriented. 

For a tetragonal lattice of symmetry classes J^mm, 422, 42m, J^/mmm, the matrix 
[R] has the form [TO] 



[R] 



( 


r\\ 


r\i 


^13 








\ 




''12 


r\\ 


»~13 















TVi, 


^33 






















''44 






















'44 





\ 

















r*66 / 



(8) 



similar to that of the tensor of elastic stiffness constants, [c]. The difference between 
the trigonal lattice of symmetry 32 and the chosen classes of tetragonal lattice is that 
ri4 and cu are zero. In that sense the hexagonal and tetragonal lattices are similar 
to each other. However, in contrast to both hexagonal and trigonal 32 systems, for 
the tetragonal lattice the condition r§Q = (r*n — ry^jl does not hold anymore. As 
a consequence, the number of independent constants is equal to that for trigonal 32 
lattice, whereas the structure of [R] is similar to that for the hexagonal lattice. 

For concreteness, we shall consider the class 422, with a system of coordinates 
chosen such that the z and x axes are the 4-fold and 2-fold rotational symmetry 
axes, respectively, while the y axis is perpendicular to both x and z. Solving the 
Christoffel equation we find that the crystal can sustain pure transversal waves in 
all three directions, x, y, and z. For the transversal waves propagating in the x 
direction, the ones polarized in the z direction propagate with the velocity yj C44/ p 
whereas the ones polarized in the y direction propagate with the velocity yj cqq/ p. The 
transversal waves propagating in the y direction are similar to the ones propagating in 
the x direction: the waves polarized in the x direction have a sound velocity of y/ cqq/ p, 
whereas the ones polarized in the z direction have a sound velocity of yj C44/ p. Finally, 
the transversal waves propagating in the z direction have all the same sound velocity, 
\AWp- 

Averaging over the directions t we obtain for the transversely polarized waves 
propagating in the x,y and z directions the results 



(|A4*,y, 



= (\M k 



ky,x,t\ 
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and 

(I 



ky,z,t\ 



= (\M, 



fez.x.t I 



2 ) = <|M fc 



a,y,tl 



4N 2 k 2 rj 4 
15 



(9a) 



(9b) 



in obvious notations: the first subscript indicates the propagation direction and the 
second denotes the direction of polarization. Plugging equations (|9"a|) and (|9bj) into 
(|7| we obtain the general type of expression 

-1 



tanh 



2k B T 



(10) 
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where the superscript *- T ) comes from the tetragonal symmetry. Concretely, the new 
7 constants are 

(7gke) 2 = {lg\, t ? = , (Ha) 

(7g,, t ) 2 = (7gk t ) 2 = (7£, t ) 2 - = (Hb) 

The ratio r§§jr^ can be obtained by measuring the characteristics (sound 
velocity change or attenuation rate) of elastic waves propagating along the (100) 
crystallographic directions. It is interesting to note now the similarity between 
relations like t^X Jt^X^ = {l^yjl^t? = (WR4) 2 and relations like 
(cfcx,y,i/cfcx,z,t) 2 = C66/C44 that hold for tetragonal lattices-e.g. Cfc x ,y,t is the velocity 
of a transversal sound wave propagating along the x direction and polarized along the 
y direction. 



3. Conclusions 



We calculated the average scattering rates of pure transversal phonon modes on the 
TLSs in tetragonal disordered crystals, to emphasize the anisotropy of the glass-like 
properties imposed by the lattice anisotropy. Two of the parameters of the model- 
the TLS-strain field coupling constants-may be obtained by measuring attenuation 
rates or sound velocity changes of elastic waves propagating along the different (100) 
directions. Determining the coupling constants enables one to calculate the glassy 
properties of the crystal in any direction. 

The calculations can be extended easily to disordered crystals of any symmetry. 
Moreover, although we used in our calculations an isotropic distribution over the TLS 
orientations, the comparison of our calculations with experimental data would enable 
one to find if our assumption is true or not. If it is not true, one can determine, at 
least in principle, the distribution of the orientations of the TLSs. 
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